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ABSTRACT 

We report on the detailed radio status of the M 87 jet during the Very-High- 
Energy (VHE) 7-ray flaring event in April 2010, obtained from high-resolution, 
multi-frequency, phase-referencing VLBA observations. We especially focus on 
the properties for the jet base (the radio core) and the peculiar knot HST-1, 
which are currently favored as the 7-ray emitting sites. During the VHE flaring 
event, the HST-1 region remains stable in terms of its structure and flux density 
in the optically thin regime above 2 GHz, being consistent with no signs of 
enhanced activities reported at X-ray for this feature. The radio core shows an 
inverted spectrum at least up to 43 GHz during this event. Astrometry of the 
core position, which is specified as ~20 R s from the central engine in our previous 
study, shows that the core position is stable on a level of 4 R s . The core at 43 
and 22 GHz tends to show slightly (~10%) higher flux level near the date of 
the VHE flux peak compared with the epochs before/after the event. The size 
of the 43-GHz core is estimated to be ~17 i? s , which is close to the size of the 
emitting region suggested from the observed time scale of rapid variability at 
VHE. These results tend to favor the scenario that the VHE 7-ray flare in 2010 
April is associated with the radio core. 

Subject headings: galaxies: active — galaxies: individual (M 87) — galaxies: jets 
— gamma rays: galaxies — radio continuum: galaxies 
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Introduction 



The location and the physical properties of Very- High- Energy (VHE) 7-ray emission 
from relativistic jets are one of the most intriguing questions in astrophysics. The nearby 
radio galaxy M 87 is a well-known VH E 7-ray emitter since the first identification of 
TeV emission by HEGRA in 1998/1 999 (lAharonian et al.l 120031 ). Thanks to its proxim- 



ity (D = 16.7 Mp c, z = 0.00436: 1 Jordan et al.l l2005h and a large bl ack hole mass (Mbh 



(3 - 6) x 10 9 M ; iMacchetto et alJll997t bebhardt fc Thomaslbood lH. the jet s tructure 



can 



be resolved under 100 Schwarzschild radii scale (R R ) with VLB! observations (jJunor et al. 



199ft iLv et al.ll2007t lAsada fc Nakamural 12012 ; lHada et al.ll201lL hereafter Hll), providing 



a unique opportunity to probe the connection of the VHE 7-ray with the relativistic jet by 
isolating the detailed substructures. 

Recently there have been three remarkable VHE flares from M 87; the events in 2005, 
2008 and 2010. In the 2005 event (lAharonian et al.l 120061 ). the VHE flare was accompa- 
nied by radio-to- X-ray flares from HST-1, a peculi ar knot located at a de-projected distance 
of at least ~120 pc downstream of t he nucleus (Harris et al. 2006 1^1. with the emergence 
of superluminal (~4c) radio features (jCheung et al.ll2007|). These lead to the strong argu- 



ment that the VHE emi s sion originates in HST-1 (e.g., IStawarz et al.l 120061 ; ICheung et al. 



20071 ; lHarris et al.l 120081. 120091). althou gh there are still some debates on this interpreta- 



tion (e.g., iGeorganopoulos et al.l 120051 ). In the case of the 2008 event, on the other hand, 



the Chandra X-ray Observatory detected an enhanced X-ray flux from the nucleus, while 
HST-1 maintained a comparatively constant flux. In addition, synchronized VLB A ob- 
servations at 43 GHz revealed a strong flux increase from the radio core that lasted the 
subseque nt ^two months. T hese provide evidence that the VHE flare in 2008 originates in 
the core ( lAcciari et al.l 120091 ) . 



The third one occurred more recently, in April 2010, where the VHE flare has been 



clear l y detected during the joint monitoring cam paign by H.E.S.S, VERITAS and MAGIC (lOng fc Mariotti 
20101 ; lAliu et al.ll2012t lAbramowski et al.ll2012l hereafter A12). The detected flare displays 
a smooth rise and decay in flux with a peak around MJD 55296 (2010 April 9-10), reaching 
a historic high state of about 20% of the flux of the Crab Nebula. Interestingly, Chandra 



1 In this paper we adopt Mbh = 6 x 10 9 M q along with lHada et al.l (|201ll ). although we note that the 
exact value of Mbh in M 87 is still controversial and should be carefully considered. One can re-scale the 
values in R s unit in this paper by multiplying a factor of 2 if Mbh = 3 x lO 9 Af0 is used. 

2 The exact distance of HST-1 depends on the viewing angle of t he M 87 jet. The o bserved superluminal 
motions of HST-1 suggest a smaller viewing ang le close to ^15° (Biretta et al. 19991) . The recent optical 
polarization study of HST-1 supports this trend (jPerlman et al.ll2011l ). 
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observations taken ~3 days after the peak of the V HE flare detected a n enhanced flux from 
the nucleus, whereas HST-1 remained in a low state ([Harris et al.ll201ll ). Short time scales of 
variabilities observed at VHE and X-ray suggest the size of the emitting region to be of the 
order of a few light days times the Doppler factor S, wh ich corresponds to les s than ~105 R s . 
We note, however, that a detailed study on HST-1 by iGiroletti et al.l (120121 . hereafter G12) 
confirmed a recurrence of this structure and its possible connection with the VHE activity. 

While the broad band light curve is studied in A12, in this letter we report on the 
detailed radio status of the M 87 jet during the 2010 VHE flaring event obtained from 
high-resolution VLBA observations. We especially focus on the multi-frequency observations 
which successfully synchronized with this event, providing a wealth of information for HST- 
1 and the inner jet during the flare on milliarcsecond (mas) scale. The data and the analysis 
are described in the next section. We then show the results in section 3. In the final section, 
we discuss the obtained results and make a summary. In the present paper, the spectral 
index a is defined as S(u) oc u +a . 



2. Observations and Data Reduction 

M 87 was observed with the VLBA at 2, 5, 8, 15, 22 and 43 GHz on April 8 and 18 
2010 (MJD 55294 and 55304), which are just before and after the date of the maximum 
flux of the VHE flare (MJD 55296). These data, which were not presented in A12 except 
for the results of the HST-1 flux at 2 GHz, are identical to those presented in Hll, where 
we investigated the core shift of M 87 using the phase-referencing technique relative to the 
nearby radio source M 84. The details of the observations and the data reduction processes 
including the astrometric analysis are described in Hll. For this paper we partly re-analysed 
the data in order to properly examine the radio status of M 87 including the HST-1 region, 
which is located ~900 mas away from the phase-tracking center. The data were averaged 
only to 5 s in time and kept individual channels (1 MHz width) separated before imaging 
process to avoid time/bandwidth smearing effects at the location of HST-1. 

In addition, we also analysed many available VLBA archival data at 43 and 22 GHz 
to investigate the lightcurve of the inner jet region around the VHE flare (see Section 3.2). 
These consist of the data on January 18 (22 GHz), April 4 (22 GHz), May 1 (22 GHz), May 
15 (22 GHz) and May 30 (22 and 43 GHz) 2010, which were not included yet in A12. 

Images were created in DIFMAP software with iterative phase/amplitude self-calibration. 
Several weighting schemes were used depending on the target region. In Figure [1] we sum- 
marize representative images of the M 87 jet obtained by our observations. 
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3. Results 

3.1. The HST-1 Region 

The HST-1 region was detected on both epochs at 2 and 5 GHz on a level of 12a and 7a 
image rms, respectvely. At 8 GHz, the analysis combined for both epochs with a relatively 
strong Mf-tapering detected this region with ~8cr. The feature was not detected at 15, 22 
and 43 GHz due to the image sensitivity limit. A part of the HST-1 properties in these 
epochs (i.e. the distance from the core and the position angle of the HST-1 region at 2 GHz) 
is already reported in G12 in the context of the long-term kinematic study of this feature. 

The overall structure observed in these epochs is similar to that on 2010.07, the image 
for which is presented in G12. The HST-1 region is resolved into two main subfeatures with 
an overall extension of ~40 mas at 2 GHz. Model fitting with two Gaussian components 
yields the sizes of these features as ~18 mas (1.5 pc, 2.5 xlO 3 R s ) and ~16 mas (1.3 pc, 
2.2 xlO 3 Rs) for the upstream/downstream components, respectively. While the emergence 
of a new component from HST-1 upstream is discovered in the later epochs of 2010 (G12), 
such a feature is not found yet in our observations on April 8 and 18. The brightness 
temperature for the brightest component is estimated to be ~1 x 10 7 K, which is similar to 



the upper limit of 9 x 10 6 K derived in the previous study at 15 GHz (IChang et al.l 120101 ) 



In Figure [21 we show the integrated radio spectrum of the HST-1 region in these epochs. 
The HST-1 region shows a steep spectrum with an averaged spectral index of a ~ —1.2, indi- 
cating that the emission region is optically thin. This radio spect ral index seems to be slightly 



steeper than the values for the optical bands (ao-uv ^ —0.7; iPerlman et al.l 120111 ) . which 
were measured between 2002 and 2007. The spectral shapes are quite similar between these 
two epochs and no significant flux variation was found within the errors. The magnitudes 
of the flux densities in these epochs appear to follow the long-term, monotonically-decaying 
trend that continues since the maximum phase in 2005. These observational characteristics 
at the radio freque ncies are consisten t with no signs of enhanced activities in HST-1 at X-ray 



during 2010 (A12: Harris et al.ll2011[ ). 



3.2. The Core and the Inner Jet Region 

The inner jet region was clearly detected at all frequecies during the VHE flaring phase. 
This region is characterized by the compact core with the edge-brightended structure as 
seen in the 43-GHz image in Figured] Similar images are obtained when images are created 
separately for April 8 and 18. 
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In Figure |3l we show the light curve of the inner jet region at 43 and 22 GHz around 
the 2010 flaring event. This is an updated version of the 43-GHz light curve presented in 
A12, where the data points analysed here were not included yet. Based on lAcciari et al. 



( 120091 ) and A12, the fluxes for three different regions are provided; (i) the peak flux when 
convolved with a beam of 0.43 x 0.21 mas in P.A. — 16° at 43 GHz or 0.54 x 0.27 mas in 
P. A.— 10° at 22 GHz, (ii) the nucleus (the deconvolved flux in the circular region with radius 
1.2 mas = 0.1 pc, centered on the intensity peak), (iii) the flux integrated along the jet 
between distances of 1.2 and 5.3 mas from the intensity peak. We assign 5% errors for all of 
the data points because the amplitude calibration of the VLBA is typically accurate within 
this uncertainty. As seen in Figure |3l the inner jet region shows a relatively moderate flux 
evolution around the VHE flare. Near the date of the VHE flux peak, the innermost region 
(i) (and also (ii)) tends to be at a slightly higher flux level (~10% both at 22 and 43 GHz) 
compared with the other epochs. This is in contrast to the event in 2008, where the 43- 
GHz core underwent a remarkable flux increase (up to ~30%) l asting the subseque nt ~two 
months together with a flux enhancement from the X-ray core (lAcciari et al.ll2009l ). 

In Figure HI we next show the radio spectra of the inner jet on April 8 and 18 for 
the three regions. Note that, regarding the peak fluxes at 22 and 43 GHz, the data have 
all been convolved with a 22-GHz beam to match the spatial resolutions. The spectra of 
each region looks the quite similar between April 8 and 18. We found the following spectral 
characteristics as a function of the measured region; the spectra gradually change from steep 
to flat toward the upstream side (a ~ —0.7 at (iii) to ~ —0.1 at (ii)), and the spectrum of 
the innermost region (i) becomes slightly inverted between 22 and 43 GHz with a ~ 0.1, 
indicating that region (i) is optically thick at these frequencies. This is consistent with the 
detection of the core shift reported in Hll, because this effect occurs only when the core 
represents an opaque part of synchrotron emission at each frequency. For some of other 
radio sources, the i r inverted ra dio cores are caused by foreground free-free absorpt ion (e.g., 



Kellermannl 11966c IO'dealll998l ). but that does not seem to be the case for M 87 (ILy et al. 
20071 ). The measured frequency dependence of the core shift j,-°- 94±0 - 09 by Hll is in good 



agreement with v 1 , which is typically expected when the co r e is dominated by synchrotron- 
self-absorption (e.g.. iBlandford fc Konigllll979l ; lKonigllll98ll ; lLobanovill998l ). 



In Figure we finally show the difference of the core shift (in RA direction) between 
these two epochs (i.e., r corCtU (t = MJD 55304) — r coretU (t = MJD 55294)). Because the 
detected core shifts on April 8 and 18 are quite similar each other, we did not find any 
significant time evolution of the core position over the all observed frequencies. The result 
at 43 GHz, where the measurement attains the highest position accuracy, indicates that 
the core position remains stable on a level of 4 R s (projected scale on the sky) during this 
period. A similar result is reported during the VHE flare in 2008 (not more than ~6 R s ; 
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Acciari et al.ll2009l l 



4. Discussion and Summary 

The location and the physical properties of the VHE 7-ray emission of M 87 are still 
under hot debate. The core and HST-1 are currently favored as the emitting sites based 
on the multi-band correlations detected in the 2005 and 2008 events respectively. In the 
case of the 2010 flare discussed here, Chandra detected an enhanced X-ray flux from the 



nucleus during the VHE flare, while HST-1 remains low state ([Harris et al.ll201ll ). This is 
reminiscent of the 2008 case. The observed timescales of short variability at VHE/X-ray 
suggest the size of the emitting region to be less than ~10<5 R s . These lead to the argument 
that the 2010 VHE flare probably originates in the innermost jet region, at least within the 
resolution element of Chandra (0'/6 ~ 50 pc; A12). 

At the radio bands, we did not find remarkable activities both in the inner jet and HST- 
1 during this event. However, the obtained results tend to favor the above interpretation 
rather than the HST-1 origin. For the HST-1 region, we did not find any compact sub- 
components that could account for the small emitting volume suggested from the rapid 
VHE/X-ray variabilities. From the observed optically thin radio spectra of HST-1, one 
can expect simultaneous correlation between radio and VHE if the 2010 event originates in 
this feature, but instead the radio flux is constantly decreasing. The discovery of a new 
component at the HST-1 upstream in 2010 is particulary intriguing (G12), but the timing 
of its appearance seems to be slightly offset from the VHE event. The radio core, on the 
other hand, remains very compact during this period. Gaussian model fitting to the 43-GHz 
visibility data indicates a deconvolved size of the radio core as ~0.12 mas (17 R s ), which 
is close to the size of < 105 R s . Moreover, the observed optically thick nature of the radio 
core is possible to hide or weaken the VHE-related activity at the radio bands, which could 
be related to the relatively moderate evolution of the radio core flux. A possibility that 
the jet between the core and HST-1 is the 7-ray emtting source is not completely ruled out 
because Chandra does not resolve this region. However, this scenario seems to be problematic 
because this region is relatively extended compared with the suggested emitting volume, and 
maintains the optically thin radio spectra with the constant flux density during the VHE 
flare. 

It should be noted that the core shift measurements obtained in these observations 
locate the 43-GHz radio core position as ~20 R s from the central engine (Hll). This implies 
that the site of the VHE 7-ray associated with the core is not more distant than ~ 20 i? s 
from the black hole. In this case, the ambient radiation field such as from the accretion flow 
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could abs orb the emitted TeV photon s due to the process of photon(VHE)- photon(IR) pair 
creation ( Neronov &: Aharonianl 2007 ) . However, M 87 has only a weak IR nucleus {yL v ~ 



10 



40-41 



erg s 



Perlman et al.ll200lf ). so this all ows the 7-ray photon s up to ~20 TeV to 



escape even from within 20 R s of the black hole (IBrodatzki et al.ll201ll ). 



If the VHE flare in 2010 really originates in the core, one should explain the reason 
for the distinct behaviour of the radio-to-7-ray correlation between 2008 and 2010. One 
possibility is the difference of the opacity at the radio bands between 2008 and 2010 due to 
the change of magnetic field strength. Assuming that the emitting plasma is spherical and 
uniformly magnetized, one can estimate the magnetic field strength through the condition of 
synchrotron-self-absorption (SSA) as B = 3.2 x 10~ 5 9 4 u^ n S^S(l + z)' 1 Gauss, where 9, u m , 
S m and 5 give the angular size of the e mitting region in mas, the SSA turnover frequency 
in GHz, the flux density at v m in Jy (IKellermann et al.lll98ll ). Adopting 9 ~ 0.12 mas, 
v m > 43 GHz, S m ~ 0.7 Jy, 8 > 1, B is estimated to be > 2.0 G. In the case of the 2008 
flare, on the other hand, a time-dependent modeling of the 43-GHz core light curve based o n 
SSA indicates a relatively moderate magnetic field strength as ~0.5 G (lAcciari et al.ll2009l ). 
If we assume that the other parameters maintain roughly the same values between 2010 and 
2008, this yields v m <35 GHz, resulting in the 43-GHz core to be partially optically thin in 
2008. Thus, more amount of the radio emission is possible to escape from the 7-ray emitting 
site, leading to the stronger radio/VHE correlation in the 2008 case. 

The radiative cooling due to synchrotron emission is also possible to contribute the 
weaker activity of the radio core in 2010 because of its relatively stronger magnetic field. 
The time scale of the synchrotron cooling at a frequency i/& can be e stimated as t synr h 



B 3 { 2 ^ , where t S ynch, B and i/& are measured in year, Gauss and GHz ( IScheuer fc Williams 



19681 ). The cooling time scale of the 43-GHz emission £ sy nch,43GHz under B > 2.0 Gauss results 
in t sync h < 20 days. This time scale is comparable to the duration of the possible decaying 
pattern of the radio core flux (seen for the region (i) between MJD 55295 and MJD 55332). 

While va rious 7-ray production models for the black ho l e vicinity/jet formation region 



2008; Barkov et al. 



are p r oposed (|Reimer et al.ll2004l ; iNeronov &: Aharonianll2007l ; iLenain et al.ll2008l ; iTavecchio &: Ghisellini 



20101 ). it is not easy to discriminate such models in these observations 



alone. To constrain the exact location and match the 7-ray e mission process to a specific 
model, the use of higher frequency VLBI (IDoeleman et al.ll2012l ) is promising, which provides 
higher transparency to the 7-ray emitting region with an event horizon scale resolution. 
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Fig. 1. — A summary of VLBA images of the M 87 jet during the VHE 7-ray flare in April 
2010. The main (global) image at 2 GHz was obtained by combining the data on April 8 
and 18 2010. The bottom right inset indicates a close up view toward the HST-1 region. 
The nomenclatures of two main features (compl and comp2) are based on Gil. The upper 
left inset indicates a 43-GHz image for the core and the inner jet (also obtained by averaging 
both epochs). The beam sizes at 2/43 GHz are 7.5x3.9 mas in P.A. —5° (bottom right in 
the inset of HST-1) and 0.43x0.21 mas in P.A. —16° (bottom right in the 43-GHz image), 
respectively. For each image, contours start from —1, 1, 2... times 3a image rms levels 
(3a = 1.0/3.3 mJy beam -1 at 2/43 GHz) and increasing by factors of 1.4. 
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Fig. 2. — Integrated spectra for the HST-1 region. The blue/red circles at 2 and 5 GHz 
indicate the data on April 8/18 2010, respectively. The black circles indicate the flux density 
at 8 GHz and the upper limit at 15 GHz, which are obtained from the averaged image for 
the two epochs. A solid line indicates the averaged spectral shape for the two epochs. For 
reference, VLBI integrated flux densities in earlier epochs (grey circles) are provided from 
the literatures: 2005.82 and 2009.64 at 1.7 GHz, 2010.24 at 5 GHz, and 2005.85 at 15 GHz. 
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Fig. 3. — 43- and 22-GHz lightcurves of the central region of M 87 near the VHE flare in 
April 2010. The fluxes for three different regions (i), (ii) and (iii) are provided (see text 
for the detailed descriptions). The light red area indicates the period where the VHE 7-ray 
flaring event was covered by the VHE telescopes (see the top panel, which is taken from 
A12). The dark red area indicates the date of the VHE 7-ray flux peak (MJD 55296). For 
reference, we show the ranges of the flux densities at 43 GHz before the 200 8 flare for the 



regio ns (i) and (ii), which represent typical flux ranges for the quiescent phase ( lAcciari et al. 
2009h 
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Fig. 4. — Radio spectra of the central region of M 87 on April 8 and 18 2010. The peak fluxes 
at 22 and 43 GHz were estimated with the synthesized beam at 22 GHz (0.54 x 0.27 mas in 
P.A. -10°). 
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Fig. 5. — Astromet ry of the core pos itions, (top panel) Detected core shifts in RA direction 
on April 8 and 18 (IHada et al.ll201ll ). (bottom panel) Difference of the core position (in RA 
direction) between April 8 and April 18 at each frequency (the positions on April 18 minus 
the ones on April 8). Positive direction indicates the jet direction. 



